Strictly correlated uniform electron droplets 
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We study the energetic properties of finite but internally homogeneous D-dimensional electron 
droplets in the strict-correlation limit. The indirect Coulomb interaction is found to increase as a 
function of the electron number, approaching the tighter forms of the Lieb-Oxford bound recently 
proposed by Rasanen et al. [Phys. Rev. Lett. 102, 206406 (2009)]. The bound is satisfied in three-, 
two-, and one- dimensional droplets, and in the latter case it is reached exactly - regardless of the 
type of interaction considered. Our results provide useful reference data for delocalized strongly 
correlated systems, and they can be used in the development and testing of exchange-correlation 
density functionals in the framework of density-functional theory. 
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I. INTRODUCTION 

Strongly correlated materials have attracted tremen- 
dous interest across different fields of physics.i- Fa- 
mous examples of strongly correlated systems are high- 
temperature superconductors, organic conductors, ultra- 
cold atoms, and semiconductor quantum dots. These 
systems provide a particular challenge to theorists - sim- 
ply because their properties cannot be predicted from the 
behavior of individual particles. 

A physically important quantity in a quantum system 
is the magnitude of the indirect particle-particle inter- 
action. This corresponds to the energy difference be- 
tween the expectation value of the quantum mechani- 
cal interaction operator and the classical interaction en- 
ergy of charged particles [see Eq. ^ below]. Lieb^ 
showed that this quantity has a rigorous lower bound 
for Coulomb- interacting three-dimensional (3D) systems. 
Later on, the bound was tighteneci^i^ and extended to 
two-dimensional^ (2D) systems. In density functional 
theory (DFT), this bound has been extensively used for 
building and testing approximations for the exchange- 
correlation energy functional (see, e.g., Refs.0-[land ref- 
erences therein). 

More recently, using physical rather than formal ar- 
guments, an even tighter bound for 3D and 2D systems 
has been proposed together with an extension to one- 
dimensional (ID) systemsi^ The basic idea of Ref. d is 
that the tightest form of the lower bound on the indi- 
rect interaction in D dimensions should correspond to 
the amount of correlation in the infinite D-dimensional 
homogeneous electron gas (HEG) in the low-density 
limiti^i^ This physically appealing idea provides an im- 
proved bound for 3D systems, the introduction of a rela- 
tively tighter bound in 2D, and a proposal for the bound 
in ID. 

Odashima and Capelle^*^ have shown through exten- 
sive numerical studies that finite electronic systems are 
energetically far above the lower bound, even when con- 
sidering the tighter form of Ref. This has triggered 



our interest to construct a finite, yet physically simple 
system that is as close as possible to the bound of Ref. 
or, if possible, even below (which would imply violation 
of the proposed lower bound). To challenge the bound 
maximally for a given density, the strict-interaction limit 
of DFT provides a suitable methodology. The mathe- 
matical structure of this approach - corresponding to a 
system with a given density and maximum spatial corre- 
lation between the electrons - has been uncovered in the 
last three yearsi^"— Consequently, explicit solutions, at 
least for centrally symmetric densities, have started to 
become availablcj^^^ 

In this paper we use the strong-interaction limit of 
DFT to investigate the Lieb-Oxford bound in 3D, 2D, 
and ID. We take the simplest imaginable test system, 
i.e., a finite Z3-dimensional electron droplet of a uniform 
density (up to a certain radius above which the density is 
rigorously zero) and examine its properties as the number 
of particles changes. Our analytic and numerical results 
show no evidence of violation of the lower bounds pro- 
posed in Ref. but in all dimensions the low-density 
result of the HEG is approached as a function of the 
electron number N. In ID, our large- A'^ limit exactly cor- 
responds to the proposed lower bound - regardless of the 
type of electron-electron interaction examined (contact, 
soft-Coulomb, and regularized). In 2D and 3D, on the 
other hand, the large-TV result is ~ 2% off the bound, 
although this small difference is within the errors associ- 
ated with our numerical procedure. 



II. THEORY 

A. Lower bound on the indirect interaction 

We consider a system of interacting electrons described 
by the Hamiltonian 
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where T is the kinetic-energy operator, t4e is the 
electron-electron (e-e) interaction, and t^xt accounts for 
an external local one-body potential. We can define the 
indirect (quantum mechanical) part of the e-e interaction 



as 



where 



U[n] ^IJdr' J drn{r)n{r')Vee{\r-r'\) 



(2) 



(3) 



is the classical (Hartree) interaction calculated from 
the (charge) density n{r). In Eq. ([2]), ^' = 
\['(ri(Ti, . . . , rATfTTv) is an arbitrary A''-electron wavefunc- 
tion (where ai denote spin variables) and n,p(r) is the 
density associated with it. The indirect e-e interaction 
has an important lower bound which can be expressed as 
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where D = 3, 2, 1 is the dimension. In 3D, the bound 
originally found by Lieb^ is best known as the Lieb- 
Oxford (LO) bound^ having a prefactor Cg*-* — 1.68. 
The bound has been tightened by physical, yet nonrigor- 
ous arguments to C3 = 1.44 (Ref. Q). In 2D the existence 
of the bound was proven by Lieb, Solovej, and Yngya- 
sott^ (LSY) with = 192v^ « 481.28. In Ref.i, a 

tighter bound of C2 = 1.96 was proposed. 

The bound in Eq. (g]) was originally constructed for 
Coulomb-interacting systems, where Vee = ~ 
Tjl^^. In ID, however, this type of interaction is ill- 
defined due to the divergence at \xi — Xj I = 0. In Ref. i 
it was shown that a ID bound can be constructed by 
applying a contact interaction or a soft-Coulomb inter- 
action. In Sec. nil CI below, the ID case is studied in 
detail considering three types of the e-e interaction. 
The bound of Eq. (|4]) can be equivalently expressed 
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is the local density approximation (LDA) for the elec- 
tronic exchange energy, corresponding to the exact ex- 
change energy for the HEG. Here the prefactors are given 
by A3 = 34/ 3^-1/3 /4 and A2 = 25/27r-i/V3 (for the ID 
case see Sec. lIII C]) . In the right-hand side of Eq. ([5|), the 
values obtained for Ad in 3D, 2D, and ID are 



A. = 1.96, 



A2 = 1.84, Ai = 2. 



(7) 



They have been proposed as the tightest bounds with 
the prefactors Cd given above - hence the bar symbol 



to differentiate from the functional A£)[^']. The upper 
bounds A u correspond to the low-density limit of the D- 
dimensional HEG. The physical argumentation^ behind 
the HEG result was motivated by the finding of Lieb and 
Oxford, who showed that there is a function A3(iV) which 
provides an upper bound for all systems with particle 
number equal to N (Ref. The function A3 (TV) is 

monotonic, with A3(iV -I- 1) > A3(iV), so that the most 
general bound of Eq. ([5]) is obtained by considering N — > 
00. 

In this paper we focus on the question how the LO 
bound can be challenged. In other words, how the wave- 
function \1/ in Eq. ((S]) must be chosen such that ADi^*] 
becomes as large as possible? For any class of wavefunc- 
tions with a given fixed density n(r) f Ref. [Tsh the answer 
to this question is 



maxA£,[5'] = Ad[?i] 



(8) 



where 



[n] = min 



mm{^\Vee\'f)-U[n], (9) 

>]>— >-n 



is the indirect Coulomb energy in the strong-interaction 
limit of DFT, which can be now calculated (at least for 
centrally symmetric densities) with the theory of strictly 
correlated electrons.— This quantity was also considered 
in the original proof of the boundj^'^'^'^ In the following 
section we briefiy review how the functional Woo[?^] of 
Eq. ([9]) can be constructed for a given density n(r). 



B. Strong-interaction limit 

We may define ^'"[n] as the wavefunction that min- 
imizes (\1/"|T -|- aVeel"^") - corresponding to a system 
where the interaction is scaled - with the constraint of 
reproducing the given density n{r). The scaled indirect 
Coulomb interaction = {■^"[n]\Vee\'^"[n])-U[n] = 

(Ke) ~ ^ W satisfies a set of useful exact relations within 
DFTJ^ 

Here we consider the strong-interaction limit a — 
CX), where it is sufficient to minimize the interaction 
term alone, since (^aVee) grows faster than (T) x a^^^ 
(Refs.HHo). As anticipated in Eq. (O, wc thus compute 
explicitly, for a given density n(r), the functional 



(K^)^ min(vI/|Ke|*> 



(10) 



In the strong-interaction limit of DFT the electrons mini- 
mize their interaction energy while reproducing the given 
smooth density n(r). This a — >■ cxd limit is thus different 
from the more commonly considered Wigner crystalf^i as 
here the one-electron density is fixed a priori (and can 
be very diiTerent from the one of a Wigner-like struc- 
ture, e.g., it can be the density of a weakly correlated 
system like a neutral atoniii). As discussed in detail in 
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Refs. [lllfl^fT^ . in the a ^ oo limit the relative positions 
of the electrons become strictly correlated: the position 
ri = r of the first electron determines the positions 
of all the other electrons via iV — 1 co-motion functions 
fi(r), Ti — fi(r). Thus, the probability of finding the first 
electron in the volume element dr around the position r 
is the same as finding the i^^ electron in the volume ele- 
ment dii{r) around fi(r), so that the co- motion functions 
are linked to the density through the differential equation 
n{r)dr = n(fi(r))dfi(r), which has to be solved with the 
boundary condition that the corresponding expectation 
value of the interaction operator^ii 



N-l N 

E E 

2—1 



dr ■ 



n(r)/N 



Uv)-lj{v)y 
(11) 

is minimum. 

Following Refs. [TTIfl3 . we consider here a spherical (cir- 
cular) density in 3D and 2D, for which the fi(r) can be 
constructed as follows. Given the expected number of 
electrons between and r. 



dr' S{r')n{r'), 



(12) 



where S{r) = Anr"^ and 2nr in 3D and 2D, respectively, 
the general solution for the radial co-motion functions in 
a centrally symmetric A^-electron system can be written 
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FIG. 1: Example of a section of the vectorial co-motion func- 
tions fi (r) for a uniform droplet in two dimensions with TV = 4 
electrons. As the position of the first electron changes on 
the a;-axis from a; = to a:: = ai = N^^{1), the position 
of the second electron changes along the boldface curve from 
r = 02 = N^^{2) to r — ai, the position of the third electron 
from r = a2 to r = = N~^{S), and the one of the fourth 
electron from r = 04 = N~^{4:) to r = as. 



functions satisfy group properties such that the result- 
ing electron-electron repulsion energy is invariant under 
exchange of two or more electrons, ensuring the indistin- 
guishability of particlesj^ii^ 
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r > a2k, 
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III. UNIFORM ELECTRON DROPLETS 



A. Three dimensions 



hk+i{r) 



iV-i[iVe(r) + 2fc], r < on- 



2k 



N-^[2N ~ iVe(r) - 2fc], r > aN- 



2k, 



(14) 

where Oi ~ ^{i) and the integer k runs for odd N from 
1 to (TV - 1) /2, and for even TV from 1 to (TV - 2) /2. For 
even TV, we need an additional function 



/Ar(r)=TV-i[TV-TV,(r)]. 



(15) 



Equations (fT51) - p^ determine the distances fi{r) from 
the center of the remaining TV — 1 electrons when, say, 
electron 1 is at a distance r from the center. The rela- 
tive angles between the electrons, leading to the vectorial 
CO- motion functions fi(r), are determined by numerical 
minimization of Eq. (jlip for each r. An example of such 
calculation for a 2D electron droplet of a uniform density 
with TV = 4 electrons (see Sec. IIIIBI for more details) 
is shown in Fig. [TJ As the position of the first electron 
changes on the a;-axis from x = to a; = ai = TV^^(l), 
the position of the second electron changes along the 
boldface curve from r = 02 to r = ai, the one of the 
third electron from r = a2 to r ~ and the one of the 
fourth electron from r — to r — a^. The co-motion 



We consider a homogeneous TV-electron droplet with a 
constant density n and a radius R. The density can then 
be expressed simply as 



^ ' I 0, r > R. 



(16) 



In a physical sense, this density and the corresponding 
wavefunction must be considered as limit cases (see the 
end of this section). The expected number of electrons 
between and r is now 

Neir)^N(^^yeiR-r)+Neir-R), (17) 

where 9 is the Heaviside step function. The Hartree en- 
ergy becomes now 



U\n] 



dr' / dr 
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and the LDA exchange energy [Eq. is 
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For N = 2, we can readily test the accuracy of the LDA 
with respect to the exact exchange energy, E'^^'^^{N = 
2) = -U{N = 2)/2. We find iJ^DA/^exact ^ 0.9621. 
This ratio is expected to approach unity as TV — >■ oo. 
The co-motion functions [Eqs. (IT3l)- (fT5)) ] become 



f2k{r) 



2k 



1/3 



(20) 



f2k+i{r) = 

r {fR^ + r'f\r<Ril^2k/N)y' 
I [(2 - f ) i?3 - r3] , r > i?(l - 2fc/A^)i/3, 



and for even N we have to add the last function /Ar(r) = 
|/j3_j.3|i/3 Tj^ggg co-motion functions keep the electrons 
in different spherical shells (each one containing, in the 
quantum mechanical problem, on average one electron - 
see the example of Fig.[T]), while keeping the first deriva- 
tive of the external potential continuousj^ The expecta- 
tion value of the e-e interaction in the strong-interaction 
limit can now be calculated from 

An J ' drr^n{r)V,e{rj2{r),...fN{r),n{r);R)^ 

R' Jo 



drr^Vee{rj2{r),...fN{r),n{r);R), 



where we have used the fact that integrating between 
and R is equivalenlii to integrating N times between 
and ai, where ai = A^~^(l) = RN~^/^. The function 
r2(r) denotes all the relative angles between the electrons 
as a function of r and is calculated numerically!^ Chang- 
ing variables x = r/R leads to 



{%^{R)) = 

3A^ / dxx'Vee{xj2{x),...fN{x)M^):R)- 



(22) 



Upon coordinate scaling it is easy to see that 

Vee{x, /2(x), . . . } N {x) ,Q.{x)] R) = 

^Veeix J2ix),... f Nix), nix); R^l) 
Finally, we can write Eq. ([8]) as a function of N, 

(v;-(i? = i))-37vV5 



A3 (TV) 



-35/377-2/3 (Ar/4)4/3 



(23) 



(24) 



where, as said, (VedR ~ 1)) is calculated numerically. 
Special cases are = 1 and N ^ 2 yielding ana- 
lytic expressions. For a single electron (V^ee) is triv- 
ially zero and we find A3(iV = 1) = 4(27r/3)2/3/5 w 



TABLE I: Calculated values for A3 and A2 as a function of the 
number of electrons in uniform strictly correlated electron 
droplets in three and two dimensions, respectively. 



N 


A3 


A2 


1 


1.310 


1.414 


2 


1.498 


1.556 


3 


1.550 


1.607 


4 


1.603 


1.644 


5 


1.627 


1.666 


6 


1.657 


1.679 


7 


1.672 


1.692 


10 


1.708 


1.719 


14 


1.733 


1.736 


20 


1.761 


1.743 


30 


1.784 


1.758 




FIG. 2: (Color online) Values obtained for A3 (circles) and 
A2 (squares) as a function of the electron number N. The 
dotted lines show the estimated limit values for uniform elec- 
tron droplets when A'^ goes to infinity. The dashed lines show 
the bounds A3 and A2 for three- and two-dimensional sys- 
tems, respectively, corresponding to the low-density limit of 
the homogeneous electron gas.~ 



1.310. For A^ = 2 we get {Vee{R = 1)) = 
3 (8- 2i/3r(l/6)r(4/3)/V7F) /20 leading to A3(A^ = 
2) 1.498. Both values are lower than those given for 
A3 (TV = 1) and A3(A^ = 2) in Ref. Tt'. 

Our numerical results for larger A^ are summarized 
in Table HI The results are also plotted (as circles) in 
Fig. [5J The curve intersecting the tabulated values has 
been obtained by numerical fitting of {V^{R =1)) using 
a liquid-drop model expansion, which leads for A3 ( A^) to 
the formula 



4 /27r 



2/3 



Af(A^) = -^Kr ai + a2Ar-i/3 + a3iv-2/; 



where ai 



(25) 



-0.879717, 02 = 0.153634, and 03 
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1.6 



1.55 



1.5 



< 1.45 
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FIG. 3: (Color online) Values obtained for A3 with A'^ = 1 
(bottom) and N = 2 (up) when a soft tail of the density 
(solid line in the inset) is applied. In the limit of a sharp 
edge (dashed line in the inset), corresponding to a 00, 
the results for the proposed homogeneous density droplet in 
Eq. (|16|) are reproduced (dashed lines). 



0.123195. When N goes to infinity, the fit yields a value 
^3{N — 00) « 1.92 plotted as a dotted line in Fig. 
This value is lower than the bound proposed in Ref. l9| 
for 3D systems corresponding to A3 = 1.9555 (dashed 
line). However, the difference is rather small (~ 2%) and 
it is actually within the error associated to the fitting 
procedure. Our values for A3(A'^) as well as our fitting 
curve are always below the model for A3 (iV) proposed by 
Odashima et ali^ 

It should be noted that, per se, the density given in 
Eq. (|16p corresponds to a non-difFerentiable wavefunction 
due to the sharp edge at r = _R. Therefore, it is impor- 
tant to examine whether the results above are valid when 
considering the density as a limit case of a physical den- 
sity. A simple choice would be a density profile of the 
form of a Fermi function, i.e.. 



n{r) 



const 



(26) 



where the numerator is the normahzation constant and 
the value for a determines the sharpness of the edge. The 
limit a — > 00 corresponds to the density of the form of 
Eq. (fT6|). Figure ([3]) shows the values obtained numeri- 
cally for A3 with = 1 (bottom) and N = 2 (up) as a 
function of a. From the figure it is clear that the (ana- 
lytic) values corresponding to the original (sharp) density 
are approached as a — cx); this is particularly convincing 
with A^ = 1 where we can numerically study very large 
values for a (even several orders of magnitude larger than 
those shown in the figure). Thus, applying densities of 
the form of Eq. (|16p , can be seen as a limiting case for a 
quantum system. 



B. Two dimensions 

In 2D a homogeneous A^-electron droplet with a radius 
R is defined by the density 



n{r) 



r < R 



(27) 



0, r > i?, 

and the expected number of electrons between and r is 



\2 

j 0{R-r)+Ne{r-R). 



(28) 



The Hartree energy is more tricky to calculate than in 3D, 
since using Eq. p8)) directly leads to an elliptic integral. 
However, we can use the fact that for a 2D disk with a 
radius r and a constant density n, the Hartree potential 
at the rim of the disk is Vh = Anr?'^ The Hartree energy 
is now equal to the work required to charge the disk from 
the center (?■ = 0) to the rim {r = R), 



N 



r27r 



The LDA exchange energy is 



dr rVn (r) 



8N^ 
3^' 



E: 



LDA 



25/2 ^3/2 



(29) 



(30) 



3ttR 

Again, for the special case of A^ = 2 we can compare 
with the exact exchange energy, iJ^^^act ^ _[//2. 
We find £'^DA/£;cxact ^ exchange energy 

is exact for A^ = 2. 

In 2D the co-motion functions have the same form as 
in 3D [Eqs. (pUj) and (PT|) ] apart from the change in ex- 
ponents as 3 —7' 2 and 1/3 —>■ 1/2. The expectation value 
of the e-e interaction operator in the strong-interaction 
fimit becomes 

2N I dxxV,,{x, hix), . . . fM{x)M^):R), (31) 



and after scaling of the distances we find 

(K~(i? = l))-8A^V(37r) 



K2{N) 



_25/2^3/2/(3^) 



(32) 



For A^ = 1 we get now simply A2(A^ = 1) = \/2, 
and the two-electron case yields A2(A^ = 2) = 2 — 
37r {8 + \/2 [In 2 ln(2 - ^2) - 3 ln(2 + ^/2)\ } « 1.556. 
Results for larger A^ are given in Table U and Fig. [21 In- 
terestingly, the 2D values are higher than the 3D ones 
at small A^, but at A^ ~ 15 they go below the 3D curve. 
Again, we use a liquid-drop-model expansion to fit our 
data, leading to 

Kf{N) - (61 + 62 Af-i/2 + 63 A^"') , (33) 

with 61 = -1.0814, 62 = 0.121609, and 63 = 0.129014. 
The large- A^ fimit yields A2(A^ — ^ 00) = 1.80, which is, 
also in this case, ~ 2% lower than the 2D bound A2 = 
1.84 in Ref.i. 
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— contact 

soft-Coulomb (|j.=1) 

regularized ( p=1) 




FIG. 4: (Color online) Different e-e interaction potentials con- 
sidered for one-dimensional systems. 



C. One dimension 



As mentioned in Sec. Ill A| the Coulomb interaction 
is ill-defined in ID. Various forms for pliysically reason- 
able ID e-e interaction operators have been suggested, 
and below we focus on three of them: the contact, soft- 
Coulomb, and regularized interaction. The shapes of the 
interaction potentials are visualized in Fig. |4l 

Regardless of the type of e-e interaction, the ID ho- 
mogeneous electron droplet has a density 



n{x) 



\A<R 

0, \x\ > R, 



(34) 



and the expected number of electrons between -co and x 
(here corresponding to the cumulative distribution func- 
tion) is 

.X r 0, X < -R 

Ne{x)^ dx'n{x') = \ ^x + !^,\A<R (35) 

^ N,x>R, 

The CO- motion functions fi {x) can be found in a straight- 
forward fashion. When the first electron [i = 1) is set at 
— i?, the electron i is located at aj_i = 2(i — 1)R/N — R. 
This leads to 

f N-^[N,{x) -f z - 1], .T < N-\N + 1 - z) 
\ N-^[Ne{x)-{N x> N-\N + 



and after substituting ^ we get 



(36) 



x + 2§{i~l), x<2§{l-i) + R 



X + 2§{i - 1) - 2R, X < 2^{1 - i) + R 



The expectation value of the e-e interaction in the 



strong-interaction limit can be written as 

n(x) X ^ 



dx ■ 



N 



Veei\Mx)-f,ix)\) 



I — u 



j) 



(38) 



where we have replaced the full integral by N integrals 
between —R and ai — —R + 2R/N, so that the difference 
between two co- motion functions is always fi{x)—fj{x) = 
2{i — j)R/N. Thus, we notice that the strong-interaction 
limit in ID, in the case of a uniform density, corresponds 
to the Wigner- crystal solution due to the translational 
invariance. In 3D and 2D instead, the distances between 
fi(r) always depend on r. 

Next, let us write the Hartree energy and the LDA 
exchange energy in a general form, 

dxdx'Veei\x-x'\), (39) 



U[n] 



8R^ 



E. 



LDA 



-R. J-R 

dxn{x)ex{n), 



(40) 



where the exchange energy per electron in a ID HEG is^ 
ex{n) = -^ dqVeM){^-7^\, (41) 



Vee is the Fourier transform of the e-e interaction and 
kp = 7rn/2 is the Fermi vector in ID. In the following, 
we will use the above general expressions to compute Ai 
in Eq. ([8]) as a function of N . 



Contact interaction 



The contact (or delta) interaction is defined as 



Vee{x) = ri6{x) 



(42) 



(see the dashed line in Fig. U) and its Fourier transform 
is simply Veeiq) = T], where 77 is a dimensionless constant. 
The Hartree energy [Eq. (|39|) ] becomes 



U[n] = 



and the LDA exchange energy [Eq. (HUl) ] is 



E: 



LDA 



8R' 



(43) 



(44) 



We see that, in this case, the calculation of the con- 
strained minimization of Eq. (jlOp does not have a unique 
minimizing wavefunction. Indeed, with the contact in- 
teraction the minimum value (V^) = can be produced 
with any wavefunction that prevents the electrons to be 
at the same position while yielding the assigned density. 
The strictly-correlated wavefunction is just one of those. 
Thus, we trivially obtain Ai = 2, which is independent 
of N, and coincides with the lower bound Ai of Ref. 19|. 
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2. Soft-Coulomb interaction 
The soft-Coulomb interaction is defined as 



(45) 



where /i is the softening (or cutoff) parameter. The po- 
tential is visualized as a solid line in Fig. 01 Its Fourier 
transform is Veeisi) — 2JCo(a*'?)j where is the modi- 
fied Bessel function of the second kindi^ This expression 
leads to the LDA exchange energy of the form 



LDA 



Ar2 

Ir 



dx2(l 



(46) 



In 



(—) 



where 7 sa 0.577 is Euler's constant. The leading term 
at small tt/zTV/ (2i?) agrees with the result of Fogleri^ 

The calculation of the Hartree energy leads to a tedious 
integral but finally yields an analytic expression, 



U\n] = 



— <^ M/i? - v/4+(/i/i?)2 - H^i/R) 




1 f 

3072 \R 













is 



o 



(47) 



where we also give the series expansion for small values 
of {11 /R), which is the regime of our primary interest (see 
below). Finally, the strong-interaction limit for the e-e 
interaction in Eq. (1551) leads to 



N 



N 

1>J 



1 



(48) 



Similarly to the 3D and 2D case, we may set R = 1. 
Thus, values for Ai essentially depend on TV and /i, and 
in E^^^ and {V^{R)) also through their product Nfi. 
Hence, in the following we fix Nfj, and examine numer- 
ically the behavior of Ai as a function of N. As visu- 
alized in Fig. 151 we find that increasing values for Nfi 
lead to a decrease in Ai, whereas decreasing iV/i leads 
to an asymptotic approach of Ai toward two. This ten- 
dency is in agreement with the bound Ai = 2 in Ref. H, 
where it was assumed that the Lieb- Oxford-like bound in 
a soft-Coulombic ID system has the same general form as 
Eq. (jU upon the multiplication of the logarithmic factor 
in Eq. (iiSl) . 
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FIG. 5: (Color online) Values obtained for Ai as a function 
of the electron number N with different values for Nfi, where 
fi is the softening parameter in the soft-Coulomb interaction. 



3. Regularized interaction 

As the third alternative for the e-e interaction in ID we 
consider the regularized form of the Coulomb interaction 
in ID. In particular, we take the representation of the 
Yukawa interaction in an infinite cylindrical wire of ra- 
dius 13 (Ref. [21I ) . The system is then studied in the limit 
where the (finite) range of this interaction is larger than 
any other length scale except the length of the wire. The 
resulting interaction potential in the momentum space 

:„21 



Veeiq) 



(49) 



where Ei(z) is the exponential-integral functioni^ In real 
space the interaction can be written as 

V;e(.) = ||e^W)erfc(M), (50) 

where erfc(z) = 1 — erf(z) is the complementary error 
function. The potential is visualized as a dotted line in 
Fig.m 

As can be expected, all terms required to calculate 
Ai become now rather cumbersome. The LDA exchange 
energy is 



E: 



LDA 



2 /•! 



— dqe'^" Ei(-qV)(l-<z), (51) 



where h = ivfiN / {2R). The Hartree integral is more 
straightforward to calculate in Fourier space. This leads 
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100 



This corresponds to the physically justified small-/? limit 
of the infinite cyhnder. Second, {V^{R)) is unstable for 
large R\i — j\/ {(iN), so we again use the series expansion. 

In Fig. ini we show the behavior of Ai as a function of 
N with different (small) values of Nf3. As in the soft- 
Coulombic case, decrease in the "cutoff" parameter in 
the e-e interaction leads to a tendency toward Ai = 2, 
although in this case the approach is very slow as a func- 
tion of N p. Nevertheless, the results are in line with 
other types of e-e interaction. It should be noted that for 
this regularized interaction a Lieb-Oxford-like bound has 
not been constructed or even proposed before. Therefore, 
on the basis of our results above we may suggest Ai = 2 
for the bound, thus agreeing with ID systems interacting 
through contact or soft-Coulomb interaction. 



FIG. 6: (Color online) Values obtained for Ai as a function 
of the electron number A'' with different values for A/3, where 
(5 is the cutoff parameter in the regularized electron-electron 
interaction. 



to 



U[r. 



1 f°° 
— J dqh^{q)V,e{q) 

N sinV?)f. . ^ 



TT' 



3/2 



2 



r-r' 



- F (1 1-2 ^-^ 

3/32 1,''''^'2' /32 



(52) 



where erfi is the imaginary error function and pFq is 
the generalized hypergeometric function.— The strong- 
interaction limit of the e-e interaction becomes 



N 



(3N 



erfc 



R\i-J\ 
I3N 



(53) 

Now, to calculate Ai for the regularized interaction 
using the quantities above, we have to restrict the pa- 
rameter range. First, E^;^^ in Eq. ([51]) is unstable for 
small b, and U[n] in Eq. (|5^ is unstable for small /3/R. 
Therefore, in both cases we numerically compute the se- 
ries expansions up to the second order of these quantities. 



IV. CONCLUSIONS AND PERSPECTIVES 

To summarize, we have examined the properties of 
strictly correlated electron droplets having a locally uni- 
form density. In particular, we have used the theory of 
strictly correlated electrons to test the validity of the 
lower bounds proposed in Ref. on the indirect electron- 
electron interaction in £)-dimensional quantum systems. 
We have found that the bound is satisfied in all dimen- 
sions, although it is approached as a function of the elec- 
tron number. In one-dimensional droplets the bound is 
reached asymptotically regardless of the type of electron- 
electron interaction considered in this work. In two and 
three dimensions, we obtain values being a few percent 
above the lower bounds. 

Our results can be taken as useful reference data 
for future investigations of strongly correlated systems 
in general, as well as for the development and testing 
of exchange-correlation functionals in the framework of 
density-functional theory. In future work we plan to 
investigate different densities, trying to challenge the 
bound maximally, as well as trying to provide reference 
data to construct A^-dependent boundsJ^ 



Acknowledgments 

This work was supported by the Academy of Finland 
and by the Netherlands Organization for Scientific Re- 
search (NWO) through a Vidi grant. 



* Electronic address: |esa.rasanen@jyu.fi| 
1 J. Quintanilla and C. Hooley, Physics World 22, 32 (2009). 
^ E. H. Lieb, Phys. Lett. 70A, 444 (1979). 
^ E. H. Lieb and S. Oxford, Int. J. Quantum. Chem. 19, 427 
(1981). 

G. K.-L. Chan and N. C. Handy, Phys. Rev. A 59, 3075 



(1999). 

^ E. H. Lieb, J. P. Solovej, and J. Yngvason, Phys. Rev. B 

51, 10646 (1995). 
^ M. Levy and J. P. Perdew, Phys. Rev. B 48, 11638 (1993). 

M. M. Odashima and K. Capelle, Phys. Rev. A 79, 062515 

(2009). 



9 



J. P. Perdew, in Electronic Structure of Solids '91, edited 
by P. Ziesche and H. Eschrig (Akademie Verlag, Berlin, 
1991). 

® E. Rasanen, S. Pittalis, K. Capelle, and C. R. Proetto, 

Phys. Rev. Lett. 102, 206406 (2009). 
1° M. M. Odashima and K. Capelle, J. Chem. Phys. 127, 
054106 (2007). 

M. Seidl, P. Gori-Giorgi, and A. Savin, Phys. Rev. A 75, 
042511 (2007). 

^2 P. Gori-Giorgi, G. Vignale, and M. Seidl, J. Chem. Theory 
Comput. 5, 743 (2009). 

P. Gori-Giorgi, M. Seidl, and G. Vignale, Phys. Rev. Lett. 
103, 166402 (2009). 

P. Gori-Giorgi and M. Seidl, Phys. Chem. Chem. Phys. 12, 
14405 (2010). 

P. Gori-Giorgi, M. Seidl, and A. Savin, Phys. Chem. Chem. 
Phys. 10, 3440 (2008). 
^® M. M. Odashima and K. Capelle, Int. J. Quantum Chem. 



108, 2428 (2008). 

M. M. Odashima, K. Capelle, and S. B. Trickey, J. Chem. 
Theory Comput. 5, 798 (2009). 

M. Levy, Proc. Natl. Acad. Sci. U.S.A. 76, 6062 (1979). 
^'^ M. Levy and J. P. Perdew, Phys. Rev. A 32, 2010 (1985). 
^° M. Seidl, Phys. Rev. A 60, 4387 (1999). 

G. F. Giuliani and G. Vignale, Quantum Theory of the 

Electron Liquid (Cambridge University Press, New York, 

2005). 

M. Abramowitz and L Stegun, Handbook of Mathematical 
Functions (Dover Publications Inc., New York, 1965). 
M. M. Fogler, Phys. Rev. Lett. 94, 056405 (2005). 
See, e.g., pp. 67 in M. Seidl, Habilitation thesis (University 
of Regensburg); Eq. (2.26) in Berkeley Physics Course, 
(McGraw-Hill, New York, 1965); pp. 41 in Ref. [2l|; Eq. 
(13) in M. Seidl, J. P. Perdew, and M. Levy, Phys. Rev. A 
59, 51 (1999). 



